S— wave pion nucleon scattering lengths from ttN, pionic hydrogen and deuteron data 
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The isoscalar and isovector scattering lengths (bo, bi) are determined using a unitarized coupled 
channel approach based on chiral Lagrangians. Using experimental values of pionic hydrogen and 
deuterium as well as low energy ttN scattering data, the free parameters of the model are calculated. 
Isospin violation is incorporated to a certain extent by working with physical particle masses. For 
the deuterium scattering length a n - d new significant corrections concerning real and imaginary parts 
are evaluated, putting new constraints from n~d scattering on the values of (bo,bi). In particular, 
dispersion corrections, the influence of the A(1232) resonance, crossed terms and multiple scattering 
in a Faddeev approach are considered. 
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The precise values of the isoscalar and isovector nN scattering lengths are one of the important issues in hadronic 
physics. Together with low energy ttN scattering data they determine parameters of the chiral Lagrangian which 
£f*l , allows to make predictions even below ttN threshold using chiral perturbation theory. The experimental data from 
' where (bo, b\) are usually extracted are the shift and width of pionic hydrogen and deuterium atoms. From the recent 
^vO \ measurements at PSI one deduces the elastic ir~p plus the n~p — * n°n transition scattering lengths 0,0,0,0]. Using 
in addition the measured a w -d amplitude 0,0 in order to determine the isospin even and odd combinations 
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a+ = - {a„- p + a n - n ) 
1 / 



a- = r [ a TT- P - a-K-n) , (!) 

or, correspondingly, the isoscalar and isovector scattering lengths 6o = <z+, b± = — a_, requires, however, a non- 
trivial work on the Tt~d system. This is because the impulse approximation (IA) vanishes in the limit bo = and 
the extraction of a„- n from a„-d calls for a multiple scattering treatment with the double scattering as the leading 
contribution. Also, higher order corrections as absorption and dispersion have an important effect, as has been 
extensively discussed for instance in ref. 0. It is in particular the ir~d scattering length that narrows down the value 
of (bo, bi). Although the error of a^- d is dominated by large theoretical uncertainties, the corrections on a n - d directly 
affect the values of (bo, b%). 

The determination of the pion deuteron scattering length from the elementary ones is one of the problems which 
has attracted much attention 0, 0, 0, m the past, but has also stimulated more recent studies 0, ^3 El 
IhH ITU Il6|] . In ref. @, a-^-d 1S calculated in a Faddeev approach incorporating several processes in the multiple 
scattering series as nucleon-nucleon correlations, absorption and the corresponding dispersion. It delivers, together 
with <b|, a very complete description of multiple scattering in the deuteron. Here, instead, we use the fixed center 
approximation (FCA) to the Faddeev equations. Other contributions, like absorption and dispersion correction, are 
evaluated separately, and fully dynamically, in a Feynman diagrammatic approach. This is feasible because the 
multiple scattering series is rapidly converging since the scattering lengths are small compared to the deuteron radius. 

In ref. the values of (60,61) have been calculated from the pion deuteron scattering length up to NNLO in 
chiral perturbation theory including the and C^ 1 ttN Lagrangians. Using realistic deuteron wave functions, and 
other modifications, the authors reobtain the double and triple scattering formulas in the isospin limit. 
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Typical results for the isoscalar and isovector scattering lengths, obtained recently in ref. |6j and ref. yj], are: 

(6 ,6i) = (-12 ±2 stat. ±8syst.,-895±3stat. ± 13 syst. ) -lCT 4 m~l @ 

(60,61)= (-34 ±7, -918 ± 13) -10" 4 m^ (2) 

Here, and in the subsequent results, the unit of inverse pion mass refers to the charged pion. The results of eq. 
are not in agreement with each other. The problem with a + is that it becomes a small quantity from a cancellation 
of terms of the order of a_, hence a + is difficult to determine, and the discrepancy between the two results in eq. 

indicates that the uncertainties in a + are larger than shown in eq. (J2J). Actually, in ref. 01 larger uncertainties 
are advocated from isospin violation, since the analysis is made by assuming isospin symmetry. The present study 
is formulated in the particle base and thus, isospin breaking effects from different physical masses are incorporated. 
This provides a part of the isospin violation [Tsl . flU I20L I2H 122 . |23| that has already been observed in a similar context 
in KN scattering [24| . 

Our first purpose in the present study is to carry out further calculations in the problem of ir~d scattering at 
threshold incorporating novel terms. We start with the absorption of the 7r~ in the deuteron and the dispersion 
tied to it. The latter contributes to the real part of a„- d , and in the literature a quite large correction originates 
from this source. Since high precision deuteron wave functions are at hand nowadays, and the analysis is carried 
out fully dynamically, a revision of the results from [3, la, Ui] is appropriate. After calculating the effects of the 
A(1232) excitation in the dispersion, other contributions as crossed terms are considered. Together with corrections 
of different nature from the literature, a final correction to a 7r - d is given. This enables us to parametrize the pion 
deuteron scattering length in terms of the elementary s-wave ttN scattering lengths a n N via the use of the Faddeev 
equations. We also test the model dependence of the results on the deuteron wave functions. 

The second purpose is then the application of the unitary coupled channel model from ref. [2j| to ttN scattering at 
low energies. The model provides the nN scattering lengths for the Faddeev equations for a^- d and the nN scattering 
amplitudes at low energies and threshold. First, it is tested if the model can explain threshold data and low energy 
nN scattering consistently. Then, a precise parametrization of the nN amplitude at low energies, including threshold, 
is achieved. 

This can be used to extrapolate to the negative energies felt by pionic atoms, the study of which has been one of 
the stimulating factors in performing the present work. 



II. SUMMARY OF THE MODEL FOR nN INTERACTION 



We follow here the approach of ref. |25j where the N/D method adapted to the chiral context of [2|| is applied. 
Developed for the case of meson meson interactions, the method of j2|| was extended to the meson baryon interaction 
in |23,|2^|, and ref. [2{| follows closely the formalism of these latter works. In the CM energy region of interest, from 
threshold up to around 1250 MeV, pions and nucleons play the predominant role compared to the influence of the 
heavier members of the meson and baryon octet. We have carried out the SU(3) study as in ref. [25j . which involves 
the KY,, KA and rjn channels in addition to the nN ones. At low energies, these channels are far off shell and we 
have seen that the fit to the data improves only slightly at the cost of three new additional subtraction constants. 
Therefore, we restrict the coupled channel formalism to n~p, n°n in the charge zero sector, and n + p in the double 
charge sector. The scattering amplitudes are described by the Bethe-Salpeter equation 



where the kernel V is obtained from the lowest order meson baryon Lagrangian 



V l3 (Vi) = -Ca ^ - Mj) 
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where Qj are the SU(3) coefficients evaluated in ref. 2jj, and G is the loop function of the pion nucleon propagator, 
which in dimensional regularization reads: 



. . , mf Mf — m3 + s Mf 

a(ji) + log -i + 1 1 log -4- 

fi z 2s mi 



log (a - (Mf - mf) + 2^ Q l (^)) + log (s + (Mf - mf) + 2^~s Q z 
log (-a + (Mf - mf) + 2^7s Q t (^)) - log (-a - (Mf - mf) + 2^7s Q, {y/a)) 
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where Qi (y/s) is the on shell center of mass momentum of the i-th pion nucleon system and Mj(mj) are the nucleon 
(pion) masses. The parameter \x sets the scale of regularization (/i = 1200 MeV) and the subtraction parameter a 
is fitted to the data. In order to ensure isospin conservation in the case of equal masses, the subtraction constants 
a^-p, cx n o n , and a^+ p are taken to be equal for states of the same isospin multiplct. Isospin breaking effects from 
other sources than mass differences are discussed later (by allowing the a^N to be different). 

The work of [25l | concentrated mostly in the region around the iV*(1535) resonance, which is dynamically generated 
in the scheme. The data around this region were well reproduced, although the description of the 7 = 3/2 sector 
required the introduction of the extra ttttN channel. The low energy data was somewhat overestimated in |25j although 
qualitatively reproduced. Here, however, our interest is to concentrate around threshold in order to obtain an accurate 
as possible description of the data in this region and determine, together with the pionic atom data on hydrogen and 
deuterium, the values of the isoscalar and isovector scattering lengths, with a realistic estimate of the error. 



1. Isoscalar Piece 



The chiral Lagrangian at lowest order that we use contributes only to the isovector ttN amplitude at tree level, 
but isoscalar contributions are generated from rescattering. Additional isoscalar terms emerge in the expansion in 
momenta of the chiral Lagrangian [Tsl l27l j and we take the relevant terms from into account following ref. 
|33j | . In particular, there is a term independent of q° with q the pion momentum, and one quadratic in q° , which enter 
into the potential V from eq. as 

Vij + 6ij { 71 ^ - 2C2 -Jf) 2M t ' (6) 

The on shell value of the C3 term (c^q 2 in ref. [33^ has been taken, consistently with the approach of refs. [25L |2S| 
which uses the on shell values for the vertices in the scattering equations. The Ci -combinations, in the notation of 
ref. [13, are fitted to the experiment. For a construction of a ttN potential up to higher energies, one has to regulate 
the quadratic term with C2, which we do by multiplying it by a damping factor, 

e -0 2 [(i°) 2 -™l]_ ( 7 ) 
Both cases, with the damping factor, and without are studied. 



2. The ttN — > 7T7riV channel 



The irirN channel opens up at CM energies around 1215 MeV. In the fits, we include energies higher than that, 
and therefore the 2-loop diagrams from this source should be taken into account as described in ref. j2]| . In ref. [^J , 
various functional forms for the real part of the irirN propagator have been tested, and setting it identically to zero 
resulted in good data agreement. Here we approximate it as a function constant in energy i/s and parametrize it in 
terms of the quantity 7. The imaginary part of the ttttN propagator has been calculated explicitly in ref. |25j| . At 
the energies of interest, it is small, and becomes only important at higher energies. 

The energy dependence of the ttN — * ttttN vertices has also been determined in ref. |2^| from ttN — > ttttN data. 
Fig. 12 of that reference shows that at low energies they can be well represented by constants, namely an = 2.6 m~ 3 
and 031 = 5.0 m~ 3 , which are the values we use. In the present approach, the irirN propagator with its two adjacent 
ttttN vertices provide nN — > ttN amplitudes which are added directly to the kernel of the Bethe-Salpeter equation 
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where all possible irirN intermediate states in the loop are considered. The 5V of eq. (JSJ are then added to the kernel 
Vij together with the isoscalar piece. In section llVl the constraints of 7 are discussed: The contribution from the ttttN 
channels should not exceed a small percentage of the corresponding Vij at irN threshold. 



3. Further refinements of the coupled channel approach 

Following the outline of ref. we take into account the Vector Meson Dominance (VMD) hypothesis and let the 
p meson mediate the meson baryon interaction in the t channel. This is justified by the identical coupling structure 
of the pNN coupling within VMD and the kernel V from eq. J3}, thus revealing the lowest order chiral Lagrangian 
as an effective manifestation of VMD. The p meson exchange is incorporated in the formalism via a modification of 
the coefficients CV,- in J3J — for details see [2£j • The explicit consideration of the p exchange helps to obtain a better 
energy dependence, reducin g th e strength of the amplitudes as the energy increases. 

One of the conclusions in [T3 was that the uncertainties of the (bo, b\) values should be bigger than quoted in the 
paper due to the neglect of isospin violation in the analysis. In the present work we introduce a certain amount of 
isospin violation by working in coupled channels keeping the exact masses of the particles. Although this is not the 
only origin of isospin violation [lSl l20l l2lt 122 . |23| it gives us an idea of the size of uncertainties from this source. Since 
there are threshold effects in the amplitudes, and the thresholds are different in different irN channels, this leads to 
non negligible isospin breaking effects as was shown in the case of KN interaction in [24| . 

Thus, the parameters for the fit of the s-wave amplitude in TtN scattering are the subtraction constant a from the 

(2) 

irN loop, two parameters from the isoscalar terms of the chiral Lagrangian, and 7 from the irN — > irirN loop. 
For the parametrization of the irN potential up to higher energies, a damping factor, introducing another parameter 
for the quadratic isoscalar term, is studied. In order to account for isospin breaking from other sources than mass 
splitting, different a,i for the three ttN channels will be investigated. 



III. PION DEUTERON SCATTERING 

The traditional approach to n~d scattering is the use of Faddeev equations [1,0,^3, although the fast convergence 
of the multiple scattering series makes the use of the first few terms accurate enough. On top of this there are other 
contributions coming from pion absorption, and the dispersion contribution tied to it, crossed terms and the A(I232) 
resonance, plus extra corrections which are discussed in detail in ref. 0. 



A. Faddeev approach 

We follow here the fixed centre approximation (FCA) to the Faddeev equations which was found to be very accurate 
in the study of K ~d scattering [34| by comparing it to a full Faddeev calculation [35ll3^ |. See also the recent discussion 
endorsing the validity of the static approximation in ref. |37| . The FCA accounts for the multiple scattering of the 
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TABLE I: Contributions to the multiple scattering series for Re a n - 





from [4] \m n x l 

J l TT J 


Phenom. Ham. [38] [m^ 1 ] 


OA) 


(-0.0001,-0.0885) 


(-0.0131,-0.0924) 


Impulse Approximation 


-2.14 ■ 10~ 4 


-0.02793 


Double Scattering 


-0.02527 


-0.02725 


Triple Scattering 


0.002697 


0.003489 


4- and higher scattering 


1.06 • 10~ 4 


5.4- 10 -5 


Solution Faddeev 


-0.02268 


-0.05163 



pions with the nucleons assuming these to be distributed in space according to their wave function in the deuteron. 
The Faddeev equations in the FCA are given in terms of the Faddeev partitions 

T„- d = T p +T n (9) 

where T p and T n describe the interaction of the 7r~ with the deuteron starting with a collision on a proton and a 
neutron respectively. The partitions at threshold satisfy 

Tp — tp -f~ tpGT n -\- t x GT X 
T n — - t n -f- t n GTp 

T x = t x + t° n GT x + t x GT n (10) 

Here, G is the pion propagator and t p , t n , t x the elementary s-wave scattering T-matrices of 7r _ on proton and 
neutron, ir° on the neutron, and the charge exchange ir~p <-> 7r°n, in this order. While the full Faddeev approach 
involves integrations over the pion momentum, the FCA factorizes the pion propagator to G ~ 1/r and eqns. (|10|) 
become a coupled set of algebraic equations. These equations are at the level of operators. At any place where charge 
is transferred from one nucleon to the other, the sign has to be changed due to the exchanged final state (Deuteron 
~ -l=(\pn) — \np))). Following ref. |34j we find for the ir~d amplitude density 

i / \ a P + On + {2a p a n - b|) jr - 2b 2 x a n /r 2 . , 

with ai being related to the scattering lengths et^ and the elementary U by 

l + — )a i = -—t i . (12) 

77Z jy / 47T 

The masses in eq. 1)120 have to be understood as the physical ones in each channel. 

The final ir~d scattering amplitude is then obtained by folding the amplitude density with the deuteron wave 
function as 



M d 



Ma 



dr \<p d {r)\ 2 A„- d {r). (13) 



This is the real part of a^-d that has to be modified by the corrections of the following sections. The latter will also 
provide the correct imaginary part of the pion-deuteron scattering length. 

If we keep up to the (1/r) 2 terms in eq. (| 1 1 11 and assume isospin symmetry, the resulting formula coincides with 
the triple scattering result of |17| up to 0(p 4 ) in their modified power counting. In order to show the convergence 
of the multiple scattering series in the the n~d collision we show in Table [I] the different contributions for two cases: 
First, for the experimental values from ref. Q with (bo,bi) = (—0.0001,-0.0885)™^, and second for (6q,&i) 
: (),() 1 .'! 1 . —0.0924)771"^ from the phenomenological Lagrangian of ref. [38| in eq. (0141) . We can see that in both 
cases the double scattering is very important and in the case of where bo is quite small, the double scattering is 
the leading contribution. 



B. Absorption and dispersion terms 

Pion absorption in deuterium has been studied in [39| using Feynman diagrammatic techniques. The absorption 
contribution reflects into the imaginary part of the (elastic) n~d scattering length. Its diagrammatic evaluation 
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FIG. 1: Absorption plus dispersion terms in n d scattering. 
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FIG. 2: Charge States in Absorption. 



leads at the same time to a dispersive real contribution to the ir~d scattering length. The evaluation of this latter 
contribution has been done using Faddeev approaches 0, 0, . We shall evaluate it here including extra terms from 
the A excitation, going beyond the non-relativistic treatment of the pions in 8, 9], and testing various approximations 
for the dispersive part. 

In order to evaluate the absorption and dispersion terms, a Feynman diagrammatic approach is used which offers 
flexibility to account for different mechanisms. We shall evaluate the contribution of the diagrams of Fig. ^ where 
Type B contributes only to the real part, including permutations of the scattering vertices on different nucleons and 
different time orderings as shown in Figs. [21 and [21 

On the first hand we consider the diagrams of type A and find the possibilities shown in Fig. For the purpose 
of evaluating the absorption and dispersion corrections we shall use the effective Hamiltonian |38l |40| 



H r = 4tt 



— Vdjifo + \ ^t($ x d° 
m-jr mi 



(14) 



with Ai = 0.0075, A2 = 0.053, which shows the dominance of the isovector part with A2. For the irNN vertex the 
usual Yukawa (f^NN / m -n) a ■ q t x vertex is taken. The value of A2 corresponds very closely to the final isovector 
term that we find, while the value of Ai is about twice as large. Yet, using the new values that come from our analysis 
in a first step of a selfconsistent procedure only leads to changes in the final results that are much smaller than the 
uncertainties found. 
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FIG. 3: Additional diagrams in absorption. 



The normalization of our T amplitude is such that the scattering matrix S is given by 
From the diagram A of Fig. ^we obtain in the tt~ d center of mass frame after performing the q° and q'° integrations 



(2tt) 4 J (2tt) 3 J (2tt) 3 
1111 



q 2 - m% + ie q' 2 - m 2 + ie 1° - e(l) + ie l'° ~ e(l') + it 

£fi* 2 *i' (q-q') (is) 

with q — (m^ — 1°, q) , g' = (m^ — i°, q') , 1' = —1, and e(l) refers to the nucleon kinetic energy. In eq. (|16fl . Fd is the 
deuteron wave function in momentum space including s- and d-wave (see Appendix), and the amplitude from the 
sum of diagrams of Fig. 121 is given by 



Sti t 2 t v to = 2 (4tt) 2 -L (2Ai + 3A 2 ) 2 C^L^ ~ 40.0 fm 4 (17) 
7774 V m n J 

where we have made the usual approximation that q° and q'° in the irN — > irN amplitude are taken as tn w /2 which 
is exact for Im T. 

The q • q' term in eq. (|16[) comes from the irNN p-wave vertices aq' aq = q' • q + i (q' x q) a after neglecting the 
crossed product term which does not contribute when using the s-wave part of the deuteron wave function. 

A different topological structure for the absorption terms is possible and given by the diagrams shown in Fig [3] 
The evaluation of these diagrams involves now the spin of both the nucleons 1 and 2 and one obtains the combination 

^1.^2 (18) 

which upon integration over q, q' leads to a structure of the type 

(T\i(T2jlilj (19) 

for the s-wave part of the wave function. The extra I integration, involving IJj and terms with even powers of 1 allows 
one to write 

&li(T2jklj * g 0-u(X2jl 2 Sij = - (?i<72l 2 = - l 2 (20) 

where in the last step we have used that c?i<?2 = 1 for the deuteron. The final result leads to 1/3 of the former 
contribution from the diagrams of Fig. [21 for the imaginary part. 
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The integration over the energy variable 1° in eq. I|16(l has been performed in three different ways. While for all 
three calculations the imaginary part stays the same as expected, the real part varies significantly, as we will see in 
the following. The result for the dispersive part for the diagrams of Figs. El and depends much on the treatment 
of the pion poles and the pion propagator. The choice of the wave function will have only moderate influence on the 
results. 

In a first approximation (Appl), only the nucleon pole is picked up in the ^-integration of eq. (|16|l . Furthermore, 
the energy components q° and q'° of the pion momenta in the propagators are replaced by the on-shell value of 1° 
which is m^/2. This is exact for the imaginary part of the elastic tt~ d scattering length. The imaginary part is given 
by cutting the two internal nucleon lines in the diagrams of Fig. and Fig. |3J and then putting the two nucleons 
on-shell. The pion mass in this picture is shared between the two nucleons that obtain a kinetic energy of m^/2 each 
after the absorption of the virtual pion on the second nucleon. 

In a second approach (App2), the pion poles of negative energy in the lower 1° half plane are still neglected, but 
for q° and q'° we substitute now the residue of the nucleon pole 

' ~ ( 21 ) 



(2^) -q2- TO 2 [m, - e(l)]" - q : 

This leads to new poles in the integration of eq. I|l()[) which correspond to cuts that affect one pion and one nucleon 
line of the loops in Figs. [21 and [31 From kinematical reasons the particles cannot go on-shell for these cuts. Indeed, 
if also the pion poles of negative energy are taken into account (App3), these poles cancel. 

Approach 3 (App3) makes no simplifications in the ^-integration any more, except the substitution of q° = q'° = 
m-^/2 in the elementary scattering length as in eq. (|17|l . The 9-dimensional integral of the amplitude (|16J) for the 
diagrams in Fig. |2in the formulation of approach 3 (App3) is: 



T = J d 3 l d 3 q' d 3 q F d (q + I)F d (q' + 1) qq' A Efi t 2 t V t 2 >, 

1 (2e(Z) - m w )(e(Z) + v)(e(l) -m,!^! (2e(l) - + 2u)(lo' 2 + c/(2e(Z) - m v ± lo)) 



(2tt) 

A = 



2e(l) -rriir-ie 2lolu'{lo + u')(e(l) + w)(e(Z) + w')(e(Z) ~ m * + - m * + w ') 

(22) 

It can be factorized to integrals of lower dimension by writing the sum of pion energies in the denominator of (|22H as 

oo 

W (ff) +«(«') ' ' (23) 



(I 



thus simplifying the numerical evaluation. The amplitude in eq. (|22|l . divided by 3, provides the imaginary part of 
the diagrams in Fig. |3 whereas the real part of the diagrams in Fig. |3 has a different analytical structure. The 
diagrams of Fig. |21 contribute with 36% to the real part with respect to the diagrams of Fig. |2I 



TABLE II: Real and imaginary contributions from absorption to a w - d for three different approaches. All values in 10 4 • m \ 





(Appl) 


(A PP 2) 


(App3) 


Im a^~ d , s-wave 


57.4 ±5.7 


idem 


idem 


Im a n - d , d-wave 


2.21 ±0.33 


idem 


idem 


Im a^-d s + d-wave 


59.6 ±5.3 


idem 


idem 


Im a n - d experimental 


63 ±7 


idem 


idem 


A Re a w - d , s-wave 


19.3 ±8.2 


13.6 ±8.9 


2.4 ±4.3 



Table [n] shows the result of all three approaches for the sum of the diagrams from Fig. 0and Fig. We have used 
two refined wave functions, the CD-Bonn potential in the recent version from ref. |4l| . and the Paris potential from 
ref. We take the average of the results obtained with either wave function. The difference of the results gives 

the error in Table [H] The statistical error from Monte-Carlo integrations has been kept below 0.1 • 10~ 4 • m _. The 
dispersive contribution from the d-wave has been only calculated for the CD-Bonn potential from ref. f° r the 

amplitude of approach 3, eq. I|22|l . The numerical value is 

A Re a^- dl d — wave, absorption = 0.18 • 10~ 4 • . 
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FIG. 4: Real and imaginary contribution to a 7r - d from absorption as a function of A from the monopole form factor. 

We also show the influence of the q'xq term, that stems from the c?q' c?q structure of the absorption diagrams of 
Fig. [21 It had been omitted in eq. (|16f) . since it contributes only in the d-wave — > d-wave transition. The contribution 
from this source has been calculated for the amplitude of approach 3, in eq. Ij22(l . for the CD-Bonn potential. We 
obtain: 

A<V- d , q , xq = (0.13 - i 0.38) • 10- 4 • m-l . 

We have also tested the relevant contributions of the absorption process with a Hulthen wave function in two different 
parametrizations taken from |43j | . The results for both parametrizations would lead to large errors of the order of 40% 
for the value of Im a T -j for the s-wave in Table ITT1 and even larger ones for the imaginary part from the d-wave. This 
is, because the two p-wave vertices make the absorption in the diagrams of Figs. [21 and sensitive to the derivative 
of the used wave function, and the Hulthen wave function is known to be less accurate, as has also been pointed out 
in ref. [6j. Therefore we do not use this simplified wave function in this study. 

Whereas the imaginary part in Table HT1 remains the same, the dispersive contribution from the s-wave decreases 
when going from approach 1 to 3. The fact that in (App3) it even changes sign for the Paris potential compared to 
the CD-Bonn potential is due to cancellations between terms, which by themselves are of larger magnitude. 

In all calculations a monopole form factor with cut-off A has been applied to the ttNN vertices of the absorption 
diagrams of Fig. [21 and Fig. Since A is the only free parameter involved in the calculation, we plot the dependence 
on A of the imaginary part, s- and d-wave, and the real part, s-wave in Fig. 0] The values correspond to approach 
3, eq. 

The real part from the d-wave is not plotted separately, since it is even one order of magnitude smaller than the 
imaginary part from the d-wave. Also plotted is the experimental value of Im a^-d = 0.0063 ± Q.0007m~_ taken from 
ref. Lai The imaginary part from the d-state has been amplified ten times in the figure. 

The results for the imaginary part depend very moderately on A. In Tablc[n] the values correspond to A = 1.72 GeV 
that has been used in the construction of the CD-Bonn potential We do not observe an amplification of the 

d-wave in absorption, relative to its weight in the deuteron wave function, as claimed in ref. Q . From the q • q'- 
structure of the absorption amplitude l|22|) one would expect an amplification of the d-wave, which has more weight 
at higher momenta than the s-wave (despite its small contribution to the norm of the deuteron wave function of 
4-6%). However, the correct combination of the angular momentum I = 2 of the d-wave, together with the spin of 
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TABLE III: Values of further contributions to the real part of a^- 





Diagram 


Value in 1(T 4 • m" 1 


A excitation 


Fig. [T] Type B 


6.4 


Crossed pions 


Fig. El 1 st 


-1.3 ± 0.1 (stat.) 


Crossed A excitation 


Fig. H 4 th 


9.5 ± 1.1 (stat.) 


Wave function correction (WFC), s-wave 


Fig. El 5 th 


-16.2 ± 0.1 (stat.) 


WFC, s-d interference, fi — 0, fjf — 


Fig. El 5 th 


14.4 ± 0.1 (stat.) 


WFC, s-d interference, [i = 0, p! = ±1 


Fig. El 5 th 


21.9 ± 0.1 (stat.) 



the nucleons in order to give a total spin of 1, leads to a very effective suppression of this enhancement. 

Also, the effect of rescattering in absorption has been investigated. In order to avoid double counting we consider 
all rescattering diagrams that have exactly one absorption insertion of the form of the diagrams in Figs. [3 and [21 In 
practice, this means a replacement of the two s-wave 7riV-vertices in the diagrams of Figs. |2Jand|31by a Faddeev-like 
rescattering series similar to eq. (|10|l . By doing so, there are no pion exchanges between the nucleons, which are 
unconnected to the external pions. Thus, effects, that are already contained in the deuteron wave function, are not 
double counted. The explicit evaluation of this class of diagrams results in negligible changes of the values in Table 
ILTIof the order of 1% or less. 

No interference between s and d-wave is observed for absorption and dispersion. In the next section and in the 
Appendix this issue is discussed further. 



C. Further corrections to the real part of a n - d 

The diagram of Fig. ^ where the nucleon pole is substituted by the A pole (Type B) is evaluated in a similar 
fashion as Type A. The sum of all possible charge configurations provides now 

t 2 tv t 2 , = 32 /^ 2 (4Ax-3A 2 ) 2 ^ 1Q g fm4 
9m 4 

where /^ata = 2.01 f n NN is the ttN A p-wave coupling. There is no imaginary part from this diagram as the A cannot 
go on shell. The numerical value can be found as 'A excitation' in Table ILTTI It is relatively small since the effect of 
the strong 7rAfA coupling is suppressed partly by a cancellation of the isovector part ~ A 2 from different charge states 
for the diagram. 

Another source of contribution for the real part of the tt~ d scattering length is given by the crossed pion diagram 
displayed in Fig. El first diagram. The T matrix for this process is given by 

T = I tB? Fd{q + 1} Fd(q ' + 1) E (tl h A t>2) (f? ■ q s ■ q,) 

(e + u) (e' + lj) + (e + e' + u) J + lu' 2 - m n (e' + to + w') 



2ujoj' (e' + uj) (e' + u>') (u) + uj') (m n — e — u>) (m n — e — uj') 
where the spin-isospin factor 



(25) 



E (h t 2 t[ t> 2 ) = 8(4^) 2 ^ (Xj - 3 A X A 2 ) ) ~ -6.0 fm 4 (26) 



nit 



is proportional to the isoscalar scattering length &o ~ Ai, and therefore the contribution of the diagram is very small. 
In order to show how a non- vanishing bo influences the crossed diagram, we have evaluated the amplitude in eq. H25JI . 
and the numerical result at the relatively large value of Ai = 0.0075 leads to the small contribution displayed in Table 
IIIII as 'Crossed pions'. The second and third diagram of Fig. El force the nucleons to be far off shell and are equally 
negligible. 

The crossed contribution with the A resonance in the fourth diagram of Fig. El is similar to eq. 125JI . the only 
difference being the nucleon kinetic energy e' which is substituted by Ma — Mm + Y /(2Ma), and the spin-isospin 
factor changing to 

2 



E(ti h t[ 4) = (4tt) 2 



6A 2 + H Al A 2 + | A? 





f%NA \ 




, ) 



56 fm 4 . (27) 
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FIG. 5: Crossed Diagrams. 



The diagram provides a larger contribution than the former ones, since it has a large weight in the isovector part A2 . 
This explains why the contribution of this diagram, displayed in Table ITTT1 (' Crossed A excitation'), is relatively large 
compared to the A excitation in Fig. ^ Type B, and the crossed pion diagram in Fig. [5] first diagram. 

Another type of correction is displayed in the last, fifth diagram of Fig. [3] The difference to the other corrections 
discussed so far concerns a pion that is not connected to the external pion line. Nevertheless, we are not double 
counting effects of the deuteron wave function and effects of the diagram. On the contrary, the last diagram of Fig. 
can be understood as a correction of the nucleon-nucleon interaction to double scattering, usually called wave function 
correction (WFC). The nucleon-nucleon interaction is, of course, modeled by a richer structure in terms of meson 
exchange as it is considered in the construction of the CD-Bonn potential of ref. |4l|. Yet, with the nucleons in the 
deuteron being relatively far away from each other, one pion exchange should give the right size of this correction in 
ir~d scattering. As one sees in Table ITTT1 ( 'Wave function correction (WFC), s-wave'), the contribution of the s-wave 
of this diagram is less than 1/10 of the one of double scattering, with the same sign (see Table QJ. To conclude, the 
correction induced by the fifth diagram of Fig. results in minor changes, that are of the size of 1/3 to 1/2 of triple 
scattering (see Table HJ. However, we do not include the contribution of this diagram in the determination of the 
corrections of a^-d- This is because it represents part of the non-static effects discussed in the next section. 

Besides the contributions of s and d-wave to the various corrections discussed so far, the interference between s 
and d-w&ve should be carefully analyzed. In ref. Q sizable contributions from this source were found. The spin 
structure of the two nucleons together with the angular momentum of the nucleons in the rf-state prohibits any kind 
of interference for the diagrams of Fig. [3 and |21 as an explicit calculation shows. Some explicit formulas for the 
angular structure of the interference can be found in the Appendix. 

The situation is different for the fifth diagram in Fig. [S] There, the pion acts similarly as in the one-pion-exchange, 
that mixes the small amount of d-wave to the s-wave of the deuteron wave function. From the p-wave character of 
the ttNN coupling we expect even an amplification of the interference of s and d-wave. This is indeed the case, as 
the numerical results in Table ITTT1 show. There, we distinguish between interference that leaves the third component 
of the angular momentum, \i and \j! (for incoming and outgoing state), unchanged ('WFC, s—d interference, /1 = 0, 
/i' = 0'), and the interference that involves different values of /j, and // ('WFC, s-d interference, \l = 0, // = ±1'). 
The latter implies a spin flip of one or both nucleons. There is some cancellation between the s-wave and the s-wave 
d-wave interference and the net effect is similar to what has been customarily taken as non-static effects in other 
works, as we mention in the next section. 

In addition to the diagrams considered so far we could add others of the type 
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v -— , — 



(28) 



which also contribute only to the real part. The approximate contribution A2 {m^ + in the ttN — > irN s-wave 

vertices becomes now A2 {—m^ + and this gives a factor 9 reduction with respect to the other terms plus an 

extra reduction from the intermediate nucleon propagator and we disregard them. 



D. Other Corrections 



Besides the dispersive contribution to the real part, the various crossed terms, and the A contribution, all of them 
discussed in the last section, corrections of a different nature occur for the real part of the ir~d scattering length. 
They have been summarized in ref. Q, and we follow here this work in order to discuss the incorporation of these 
effects in the present study. 

Fermi motion / Boost Correction 

The single scattering term for the irN p-wave interaction gives a contribution when the finite momentum of the 
nucleons in the deuteron is taken into account. A value of 61(7) • 10^ 4 m^ 1 arises from the study of ref. [fj tied to the Cq 
coefficient of the isoscalar p-wave ttN amplitude. For the s-wave such contributions cancel with other binding terms 
according to ref. 0]. In a different analysis in ref. a finite correction arises from the s-wave interaction which 
is tied to the C2 coefficient in the chiral expansion. The full C2 term in ref. [Sfij ] provides a momentum dependence 
which accounts for p-wave plus also effective range corrections of the s-wave. A large fraction of the ci coefficient in 
[30| is accounted for by the explicit consideration of the A in [^,[4^. The contribution to the n~d scattering length 
in ref. |17| depends much on the prescription taken for the expansion and in the NNLO expansion it gives one order 
of magnitude bigger contribution than the NNLO* expansion. The value of the NNLO* expansion is considered more 
realistic in ref. |l7j . We adopt here the scheme followed so far, looking into different mechanisms and using the Bonn 
and Paris wave functions to have an idea of uncertainties. For this purpose we evaluate the contribution due to Fermi 
motion in s and p- waves. 

By using the formalism of ref. jjfj we have up to p-waves for the ttN scattering amplitude 

T = 6 + 6l(t • r) + [c + ci(t • f)] q' • q. (29) 
In addition a range dependence of the s-wave amplitude is used in |6( where 

a v - p (u) = a+ + a- + (b+ + b~) q 2 

a„- n (u) = a + -a- + (b + - b~) q 2 (30) 

By performing a boost to the irN CM frame where cqns. H29I30|) hold, we obtain the impulse approximation for the 
p-wave contribution, including the correction for the range of the s-wave part as: 



K-d = 2 [cq + b + ) - 



m N 



(31) 



which coincides with ref. [|| up to small corrections of 0(p 4 ) and the introduction of the s-wave range parameter 
correction which amounts to a 25 % decrease of the term We obtain from eq. I|31|l a contribution of 57±9TCP 4 77i 7i r 1 
with 9 % of this value coming from the p 4 term. 

Next, we also take into account corrections due to double scattering with one s-wave and another p-wave, or two 
p-waves. This leads to small corrections but with large uncertainties which are genuine and should be taken into 
account as we show below. Another reason to explicitly evaluate these corrections is that they are included in what 
is called dispersion corrections taken in ref. from ref. 



13 



TABLE IV: Fermi corrections to a 7r - d in double scattering. 



Contribution 


Bonn from ref. [41J 


Paris from ref. [42J 


[icr 4 -™- 1 ] 


s 


d 


s 


d 


q 2 -term of eq. 1321 


-30 





-2 





Rest of eq. 


13 


9 


11 


13 


From eq. 


-3 


2 


-4 


2 


Sum 


5± 15 



We consider double scattering with one wN vertex in s-wave and the other one in p-wave and viceversa, or the two 
vertices in p-wave. After performing the boosts to the CM frames we find 



(s-p) , (l + m^/m N y 
a _ , = 47r 



1 + rmr/md 

2 



(2b c - 46 lC i) J ^3 i J d 3 p£{p) ^(p + q)* 



■in t, 



(P 2 + (P + q) S 



m N 



(32) 



with <p(p) the Fourier transform of the deuteron wave function <p(r), normalized such that J d 3 p \<p(p)\ 2 = 1. The 
large term in eq. I|32l) comes from the q 2 of the square bracket in the equation. This term is easily evaluated since it 
becomes proportional to \ip(r — 0)| 2 . The precise value of <p(0) is not well known since it depends on the short range 
forces assumed in the model. Thus one can imagine that there will be large uncertainties in this term. For the d-wave 
part of the wave function it vanishes, but the s-wave part provides a contribution, as can be seen in Tab. E| which 
depends much on the model. It is interesting to notice that in the case of the Paris potential where there is a stronger 
repulsion at short distances, the value of the correction is much smaller than that for the Bonn potential. We do not 
consider here the range of the s-wave part. It does not come as in the impulse approximation but it produces a small 
correction to another small correction. 

As for the p-wave in both vertices the contribution is 



Ap-p) _ 



47T 



(1 + m K /m N ) 



1 + inir/md 



m N 



(p + q) 



m N 



(p + q)q 



m N 



l»T, 



m N 



pq 



(33) 



This contribution should be smaller than the former one by comparing the strength of double scattering in s-wave with 
double scattering with s and p-waves. A straightforward evaluation with a monopole form factor in each vertex with 
A = 1 GeV, since beyond that momentum, the wave function is certainly unreliable, gives indeed a small contribution 
compared to the typical size of the corrections discussed. The results obtained with eq. (|33|l for the double scattering 
with p-wave, shown in Tab. IIVI agree with the value of —3 in the same units quoted in refs. |fjlll5j|. 

Isospin violation 

The topic is thoroughly investigated in 0] but not considered in f° r the evaluation of the ttN scattering lengths 
from pionic atom data. Their non consideration reverts into admittedly smaller errors in (&o, &i) than the given ones 
according to ref. ^tJ- The effects from this source are estimated relatively small in the n~d scattering length, of 
the order of 3.5 • lO^m" 1 according to 0, 0|. We shall take a different attitude. Our approach allows for 
isos pin violation since the masses of the particles are taken different. This is not the only source of isospin violation 
[l8L Il9l I20L I2H I22I |23| . but it gives the right order of magnitude. In order to consider breaking from other sources 
than mass splitting, we shall allow the subtraction constants on in the three 7rA^ channels to be different. The global 
fit from the section [Tvl takes into account the isospin breaking in the ir~d scattering length. 

Non-localities of the irN s-wave interaction 

These are corrections to the assumption of point-like interaction in the irN vertices. They are considered for single 
and double scattering in ref. [|| leading to a modification of the pion propagator G ~ 1/r. The non-locality of ttN 
interaction affects mainly the isovector part of s-wave 7rA^ scattering This is closely associated with the VMD 
assumption, which states that the 7rA^ interaction is predominantly mediated by the p — meson. In the work of ref. 
|25| the p meson is explicitly taken into account, modifying the Weinberg-Tomozawa term at intermediate energies, 
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TABLE V: Corrections to Re a„- d . 



Contribution 


Value in 10~ 4 • m" 1 

7T 


Source 


(jv~p, yn) double scattering 


-2 


mi 


Form factor / Non— locality 


23 ± 15 


\6 15 161 


Non-static 


11 ±6 


[6]([15, 16, 441) 

mm v mmm> mmm? mmm) ' 


virtual pion scattering 


-7.1 ± 1.4 


[49j 


Dispersion 


2.4 ±4.3 


Present study 


Crossed tt and A(1232) 


14.6 


Present study 


Fermi motion, I A 


57 ±9 


Present study 


Fermi m. double scatt. (s-p,p-p) 


5 ±15 


Present study 


Sum 


104 ± 24 





but close to pion threshold the p does not pla y a role in the approach of [25j and in the present work. Therefore, 
we adopt the corrections from refs. [S El where values of 17(9) • lCP 4 ?™^ 1 and 29(7) • lO^m^T 1 are obtained, 
respectively. We shall take a value of 23 • lO^m^ 1 , and a larger uncertainty of 15 in the same units, in order to 
account for the discrepancy of these two results. 

Nonstatic effects 

These are corrections that go beyond the assumptions made from fixed centers. The accuracy of the static approx- 
imation is further supported by the study of 37] where, due to the dominance of the isovector irN amplitudes, recoil 
corrections are shown to be small. Estimates of these nonstatic effects are done in 0, 0, 0] and, as quoted in , 
they lead to a correction of 11(6) • lO -4 ?™" 1 . Since the 5th diagram in Fig. [51 is part of the nonstatic effects, we do 
not include its contribution in the corrections, but adopt the value from |15L I16L l44|. 

Dispersion corrections and other real parts of the amplitude 

The dispersion corrections tied to the absorption of Figs. El and have been calculated, e. g., in refs. j^EJElEl- 
In ref. [8j, a repulsive contribution of —50 ± 3 • 10 _4 r7i 7i r 1 is found (the error taken from the precision of displayed 
decimal digits in ref. [8j ) . The authors of ref. |9j consider any absorption contribution to the real part from reactions 
of the type ir~d — > NN — > 7r~ d within their non-relativistic treatment of the pion. In the formulation of the present 
study, this would correspond to the sum of the absorption diagrams of Figs. |21 and [31 plus diagrams that contain 
external pions which couple directly via p-wave to a nucleon. The latter contribute only when Fermi motion is 
considered (see above). Their diagram D would correspond to the diagrams of Fig. [5] in the present workjand D' 
would correspond to the diagrams in Fig. [3J Using different models and wave functions, the authors of ref. Q obtain 
values for the sum of the diagrams D and D' of ±14.9, —19.5, and —9.28 in the units of 10 _4 m~ 1 . This tells us that 
there arc large intrinsic uncertainties in this calculation. 

In ref. [llj], a value for A a7I - d of -56 ± 14 • 10" " 4 m_ 1 was deduced from the two publications 0, . The value of 
ref. 0] is then adopted by the recent work of ref. . The other important diagram included in ref. @ , called C, 
which leads to this final number, can be interpreted in a Feynman diagrammatic way as double scattering involving 
one ttN vertex and a second scattering involving only the p-wave mediated by the nucleon pole. According to the 
work done here this should be complemented also with the A pole, which is dominant, and the crossed nucleon pole 
term. Altogether, this would be the Fermi correction to double scattering with s and p-waves evaluated before, and 
in some case (see Tab. II V (I we found a correction of —30 • IO^to^ 1 , which follows the trend of the result of ref. [Tl) . 
but we also discussed that this contribution is very uncertain since it depends on the unknown value of f(0) which is 
highly model dependent. 

We, hence, follow our Feynman diagrammatic technique accounting for the mechanisms implicit in refs. 0, and 
substitute their numbers by our dispersion correction plus the rescattering terms involving p-waves discussed above 
in the subsection of Fermi corrections. 

Other possible contributions related to pion interaction with the pion cloud are shown to cancel with related vertex 
contributions in ref. [l7| . something also found in K + and n nucleus scattering in refs. 0] and 0]i respectively. 

In addition, there are other minor contributions from the literature to the real part of the pion deuteron scattering 
length, and we show them in Table IVl without further comments. 
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IV. RESULTS 

The threshold data from pionic hydrogen and deuterium are obtained from the PSI experiments of 0, 0, H3, 01 ■ 
According to preliminary results from PSI experiments [50| . these results could change in the near future with a 
consequence in the extracted values of bo, b\. Systematic uncertainties were also considered in [(|. In addition, the 
Coulomb corrections on the pionic hydrogen have been recently revised in ref. [ElJ ■ This shifts the hydrogen data by 
the order of one a. We take the newer values from ref. |5l) but keep the more conservative error estimates from refs. 

0,11113 

a^-p^-p = (870 ± 2 stat. ± 10 syst. ) • lO" 4 /?^ 1 

a^-p^n = -1250(60) • lO" 4 ?^ 1 

a^-d = (-252 ±5 stat. ± 5 syst. + i 63(7)) • lO^m" 1 (34) 

The sum of the corrections to the real part of the pion-deuteron scattering length from the sections IIII Bl IIII CI and 

ITTTDl is 

(104 ± 24) • lO -4 ™," 1 . (35) 

The corrections in l|35|l are positive, which means that this additional attraction of the pion must be compensated 
for by a larger contribution from the multiple scattering series in order to give the experimental value in l|34|l . 

We proceed now to determine the isoscalar and isovector scattering length (bo,b\) from the experimental data in 
(|34J) . and from low energy ttN scattering data [56j . As Table [i] shows, the deuteron scattering length is particularly 
sensitive to bo- This is due to the impulse approximation that cancels in the limit bo — > 0, but contributes significantly 
for bo as the results for the phenomenological Hamiltonian H4|) in Table [I] demonstrate. The data from pionic 
hydrogen, a lx - v ^ lx - v and a^-p^TrOn, on the other hand, provide exact restrictions on the isovector scattering length 
b\. Finally, the low energy 7riV scattering data, together with the threshold values, determine the free parameters of 
the coupled channel approach of section ITU 

The authors of refs. [52LI53T ] use a different approach, based on effective field theory (EFT) in the framework of QCD 
and QED, which gives somewhat different values for the 7r~p scattering length, with larger errors tied basically to the 
poorly known /i parameter. That means, for instance, that these authors explicitly take into account corrections that 
go into the generation of mass splittings, while the empirical analyses that lead to the data in eq. (|34|l and scattering 
data at finite energies, that we shall use later on, only deal with electromagnetic corrections through the Coulomb 
potential and use physical masses. Hence, the simultaneous use of data extracted through these different methods 
for the global analysis that we use here, should be avoided. A reanalysis of raw ttN and ird data through the EFT 
techniques is possible, and steps in this direction are already given in [52L l53l l54l|. In any case, we shall show later 
on how our final results change with different values of this ir~p scattering length. The global analysis of the bulk of 
threshold and scattering data that we do in the present work leads us to use the phenomenological multiple scattering 
method which has been used to extract the amplitudes, which implies that the results we obtain should only be used 
within this framework. 



A. The isoscalar and isovector scattering lengths 



1. Previous Results 



In the literature, different values for {bo,b\) have been extracted from pionic atoms and low energy irN scattering 
extrapolated to threshold. In ref. [5^, a value of a+ = —80 ± 20 • lO^m^ 1 has been extracted from the data by 
performing a partial-wave analysis. Note, however, the comments in ref. [3 on this value, concerning the outdated 
database and other uncertainties. 

Extrapolations of low energy nN scattering to threshold have been updated over the years [561 l57l |5S| , and the 
value for the isospin even scattering length in the SM95 partial wave analysis is a + = —30 • 10~ 4 m~ 1 . In an earlier 
publication in ref. [5£|, the same group found a deep minimum in their global fit for a+ = -100 • lO" 4 ™^ 1 . The 
current value is a + = —10(12) • \{)~ A m~ 1 in the most recent analysis, FA02, of ref. [5q|. 

From the constraints of the strong shifts in hydrogen and deuterium pionic atoms, the authors of refs. 00 deduce 
small and positive values for a+ of to 50 • IO^to^ 1 . The findings in ref. 0| are still compatible with isospin 
symmetry, although the two bands (constraints) from the shift at a \ 2 of 1 do barely intersect with the constraint 
from the hydrogen width. As pointed out in ref. , this would be evidence of isospin violation. The value in ref. 
of (bo, bi) relies on the corrections to the real part of a^-^ from the analysis of ref. [111 ]. In a recent publication on new 
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FIG. 6: (Color online) Constraints on (bo,b\) from Re a n - d and pionic hydrogen. 

measurements of a^- d , ref. 3, one finds an extensive discussion on updated corrections, including refs. [ft ITEIIT^ lifjj . 
The authors in |3j find a value of (bo, &i) = — li|i ' lO -4 ™" 1 . See the discussion in ref. |(| on this value. The recent 
theoretical approach in ref. provides a value of around bo = (—30 ± 40) • 10^ 4 m^ 1 . 

To conclude, the experiments and subsequent analyses on pionic atoms lead to a value of bo being compatible with 
zero, or slightly negative, with errors of the same size or much larger than bo itself. In the various extrapolations of 
low energy nN scattering data to threshold, more negative values of bo are favoured. 

2. Constraints on (bo,bi) from threshold data 

The restrictions on (bo, bi) can be separately analyzed for the three threshold data points of eq. J2U, and for low 
energy irN scattering. The separation of threshold and finite energy allows for a consistency test of the data, and, 
on the other hand, for a test of the freedom that the theoretical model has, when only one threshold point is fitted. 
The influence of the data in eq. I|34|) corresponds to bands in the (bo,bi) plane, whose width is determined by the 
experimental and theoretical errors. 

We begin with the influence of a^-d on (^0i F° r the- real part of the deuterium scattering length a n - d-^ir- d the 
band in the (bo, b±) plane is calculated in four different ways, in order to determine the effect of isospin breaking by 
using physical masses instead of averaged ones. In all approaches, the large correction from eq. I|35(l is taken into 
account, and in the approaches one to three only the experimental error is considered. In the fourth approach, the 
large theoretical error from eq. (|3~5)) is added, widening the band by a factor of 3. In Fig. El the results for the 
deuterium are plotted. 

• First, the four irN scattering lengths in the solution of the Faddeev equations in eq. are expressed in terms 
of bo and b\. This, of course, implies the assumption of isospin symmetry. Then, random values of (bo,b\) are 
generated and a^-d is calculated with the help of eqns. I|11I13|) . The pairs (6o, b%), that lead to a \ 2 < 1 with 
the experimental value a^-d, cxp of eq. I|34f) . are kept, and plotted in Fig. (dotted line). 

• In a second approach, we generate the four scattering lengths with the coupled channel (CC) approach of section 
ITT1 For that, we take random values for the five free parameters of the theory (a, (3, 7, Ci), and also use averaged 
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masses for pions and nucleons, an assumption that we will drop in the third approach. Then, from the scattering 
lengths, (6 ,6i) are calculated with the help of eq. (|TJ. At the same time, a 1T - d is calculated with the help of 
the Faddeev equations I|11I13() . The same selection rule oi x 1 <1 as in the first approach sorts out the (60,61) 
pairs that are plotted in Fig. [f)]with the dashed line. The result coincides exactly with the first approach. This 
is indeed expected since we take the same subtraction constant a^N for all nN channels. 

• In a third approach, we use physical masses instead of averaged ones, and proceed otherwise exactly as in 
approach 2. The result, plotted with the thin solid line in Fig. [fjl shows a nearly identical result compared to 
the first two approaches. This gives us a measure to which extend isospin breaking effects from different masses 
can affect the result. The maximal effect of isospin breaking from this source cannot exceed a small fraction of 
the experimental error. 

• The fourth approach takes into account the large theoretical error from eq. i|35|) of 24 ■ 10^ 4 m^ 1 and follows 
otherwise the third approach. The band is widened significantly (thick solid line). In the following calculations 
we use this approach. 

The constraints from pionic hydrogen are also plotted in Fig. For each of the two bands from pionic hydrogen, 
marked as '-K~p — > 7r~p' and 'ir~p — > 7r°n' in Fig. we have followed the approaches 1 to 4 as for the deuteron. The 
four approaches give identical results for each band, as before. Therefore, in Fig. one finds only one band from 
a-n-p—tTr-p an d one from a T - p ^^o n . 

The horizontal band shows the constraint from the experimental a^-p^^On that is directly related to the hydrogen 
width. We have the isospin relation 

61 = 1/V2 ( a 7r - p _ 7r o„ ± Aa 7r - p ^ 7r o„) (36) 

where the experimental error A is relatively large, leading to a wide band. The hydrogen shift is closely related to 
a n - p ^ir-p: which leads to the constraint 

bi = bo — a n - p ^ 7T - p ± Aa^-p^n-p (37) 

in the isospin limit. 

Taking exclusively the data from pionic hydrogen, values of 60 from —70 ■ lO -4 ™^ 1 up to positive numbers are 
allowed from Fig. Then, the band from pionic deuterium is added which appears with a steep slope and narrows 
significantly the region of allowed values of 60. The range of 60 is now determined by the position and width of 
the deuterium band, namely by eqns. I(34|) an d (|35|1 . This shows the necessity of having revised and extended the 
corrections of the ~k~ d scattering length in the former sections. Indeed, if we would not have applied the corrections 
from eq. 1)35(1 . the deuterium band would show up in Fig. H3with the same slope, but shifted by around +55 • lO^m" 1 
along the 60 axis. This would lead to a value of 60 being perfectly compatible with 0. However, the situation after 
applying the corrections leaves us with a 60 G [—70, —40] ■ IO^to^ 1 . 

3. Pion Nucleon Scattering at finite energies 

The unitarized coupled channel approach is applied in order to describe irN scattering at finite energies. We fix the 
free parameters of the theory by fitting the model to the data above threshold. Then, the threshold prediction of the 
model is calculated. This is called 'Extrapolation' in the following. Comparing the predictions at threshold with the 
experimental data from pionic atoms, eq. i|34|) . the low energy behavior and the consistency of the model is tested. 

Additionally, it is desirable to have an accurate parametrization of the ttN amplitude over some energy range, from 
threshold up to moderate energies. This is achieved by including the threshold data themselves in the fit, and this is 
referred to as 'Global fit' in the following. 

Selection of experimental data: From the analyses of the CNS data base for ttN scattering we choose the 
'single-energy solutions' values, which are obtained by fitting narrow regions in the CM energy yfs separately. In 
contrast to the global fit given in ref. j5(J, the single energy bins carry individual errors each. This helps to determine 
the statistical influence of ttN scattering on the par ameters of the model. We add a small constant theoretical error of 
0.002 to the amplitudes in the normalization of [5||. The channels to be included are the s-wave isospin 7=1/2 and 
7 = 3/2 amplitudes, with real and imaginary part. This does not mean four independent data points for each point 
in energy ^/s: Since the inelasticity is zero at the low energies included in all fits, real and imaginary part are totally 
determined by the phase shift S, and therefore, for the purpose of calculating the reduced there are only two 
independent values, from the 7=1/2 and from the I = 3/2 channel. In Fig. [7| the global fit and the extrapolation 
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FIG. 7: (Color online) Threshold data from D and H, extrapolations from ttN scattering, and global fits. 



TABLE VI: Global fits and threshold extrapolations from Fig. [7| Also, the fit without the damping factor from section II V A 51 
is displayed. 





Global Fit 


Extrapolation 


No damping 


fitted data (^/s) 


1104-1253 MeV + threshold 


1104-1253 MeV 


1104-1180 MeV ± threshold 


xl 


51/(2- 10 + 3) ~2.2 


24/(2 ■ 10) = 1.2 


33/(2 ■ 6 ±3) = 2.2 




-1.143 ±0.109 


-0.990 ± 0.083 


-1.528 ±0.28 


lex - c 3 [GeV" 1 ] 


-1.539 ±0.20 


-1.000 ±0.463 


-0.788 ±0.14 


c 2 [GeV" 1 ] 


-2.657 ±0.22 


-2.245 ±0.45 


-1.670 ±0.07 


(3 [MeV~ 2 ] 


0.002741 ± 1.5 ■ 10~ 4 


0.002513 ±3.3- 10~ 4 


No f3 


7 [10- 5 • ml] 


5.53 ±7.7 


-10 ±6.1 


10 ± 10 


X (^7T — p — *7T ~ p ) 


3 


[91] 


4 


X (^7T~ p~ *7T D 71 ) 


< 1 


[2] 


< 1 


X 2 {a^- d ) 


8 


[6] 


4 



of irN data to threshold are plotted. Table IVT1 displays the parameters of the fits, including the range of energies of 
the fitted data. 

The global fit favors values for bo and b\ still negative but with smaller strength than the threshold data. As this 
fit contains threshold and finite energy data, it is situated between the extrapolation and the intersection of the three 
bands in the (60,61) plane. These differences between threshold and extrapolation from scattering data mi ght be 
related to possible additional uncertainties in the phenomenological extraction of the partial wave amplitudes jSflJ ■ 

In Fig. [3 the extrapolation and the global fit are indicated by shaded regions. These regions can be understood in 
the way that the reduced Xr from Table IVll does not raise by more than 1 from the optimum (A%^ < 1) for all points 
(bo,bi) inside these regions. For both fits, the optima are indicated with the black dots in Fig. [7] Furthermore, we 
plot for the global fit the region that fulfills A\r < 2. It appears as the light grey area just around the Axl < 1 
region. 
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TABLE VII: Isoscalar quantities: j3, isoscalar generated by rescattering in nN, and final result of the fits. 





Global fit 


Extrapolation 


b, : 


[HTW 




-336 


-434 


bo 


[kt 4 ™- 1 


, generated 


442 


396 


bo 




, final 


-28 


-46 



The above explanation for the shaded regions of Fig. 0has to be taken with caution: The values (6o, bi) inside the 
shapes have been calculated by the use of eq. Q from the elementary scattering lengths in the particle base. They 
are not the free parameters of the theory, which are the a, (3, 7, C2, and (2ci — C3). 

This implies that the A\ 2 < 1,2 criterion is applied to a x 2 that is a function of a n - pi a^- ni o T o n , and a 7T - p ^ 7T o n . 
Once a set of elementary scattering lengths a^N fulfills the criterion, (60, b\) are calculated from these values via eq. 
and give a point in the shaded regions of Fig. 

The elementary scattering lengths themselves have been calculated with the help of the CC approach by gener- 
ating the fitting parameters randomly in a wide range. Every set of a^jv corresponds to a set (a, /3,7, Cj) and, via the 
criterion, the parameter errors on (a,/3,7, Ci) in Tab. IVII are determined. 'Parameter error' means here: The range 
of a parameter of a model, that leads to a raise of the reduced \ 2 of less than 1 from the best % 2 , minimizing \ 2 at 
the same time with respect to all other parameters. 

As the final results for (60, 61) we take the values from the global fit: 

(6 , = (-28 ± 40, -881 ± 48) • [lO" 4 ™" 1 ]. (38) 

The errors have been taken from the maximal extension of the region in (60, 61), calculated from the Ax 2 < 1 criterion 
described above. They can be read off Fig. [7| The errors on (60, 61) take also into account the uncertainties from irN 
scattering data up to 1253 MeV. 

With the caveat expressed in section lTVl about using simultaneously data obtained through the EFT or phenomcno- 
logical multiple scattering methods, we would also like to give, only as indicative, the results that we would obtain if 
we replaced the scattering length of n~p from eq. I|34|l by the one given in ref. [53|. We find then the values 

(b , bi) = (-39 ± 52, -862 ± 68) • [KrV^ 1 ]. (39) 

which are compatible with the results in (|38|l within the error bars. The changes from eqns. I|38|l to (|39|) go in the 
same direction as in ref. [17| fo r bp but not for b\. This is due to the fact that our fit puts more weight in the 
scattering data than that of ref. [17j . 

4- Finite energy behaviour of the fits 

In Fig. |S| the energy behaviour of the global fit and the extrapolation from Fig. [7|and Tab. IVII is displayed. The 
two upper pictures show the real and imaginary parts of the isospin 7=1/2 channel, the lower the same for / = 3/2. 
The data from the CNS in the single-energy solutions is displayed with errors. The global fit is displayed with 
the solid line, the extrapolation with the dashed line. As expected, the extrapolation provides better high energy 
behaviour, as it is not restricted by threshold data. It is remarkable, how well the extrapolation matches the data 
above 1253 MeV, which is the upper limit of the fitted data. 

In ref. [2{| the irN scattering data has been fitted up to high energies, including the region of the N* (1535) 
resonance. A fit was obtained that explained well the resonance but overestimated the 1=1/2 amplitude at low 
energies, even when includin g th e p-meson in the i-channel (see sec. [HJ. The itttN channel does not substantially 
improve the situation in ref. |25| . It seems to be impossible to have at the same time a precise low energy fit, and a 
reproduction of the iV*(1535) resonance with the input of ref. 25]. In the present approach we have introduced the 
extra isoscalar term which improves the fit of the low energy data. On the other hand, the present approach cannot 
reproduce the N* (1535) resonance since the heavier members of the baryon and meson octet become important at 
these energies and one has to use the full SU(3) approach as in ref. |25j . 
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5. The size of the isoscalar piece and the (n, 2n) term 

In table IVlJ we compare the two sources of isoscalar strength. In the first line we show the value of the contribution 
to 60 from the isoscalar term of eq. (JBJ , 

1 m N 4ci - 2c 2 - 2c 3 2 
°c = ; 73 "V (40) 

The second line of Tab. IVHI shows the 60 that is generated by the rescattering of the irN system. It has been 
calculated for the fits by setting all parameters except the subtraction constant a to zero. In this way, one can 
extract the size of the isoscalar part that is generated by the multiple loop sum from the Bethe-Salpeter equation 
0. Although the lowest order chiral Lagrangian from ref. provides pure isovector interaction, the rescattering 
generates an isoscalar part, where usually 90 % is generated by one loop, and most of the rest by the 2-loop rescattering 
(depending on the actual values of the subtraction constants). The last line of Tab. IVIII provides the final value of 
bo- The interplay of the isoscalar piece from eq. ©, the isovector interaction, and the subtraction constant leads to 
a resulting bg (last line of Tab. IVII|I that cannot be explained any more as the sum of b c plus the generated bo- 
lt is instructive to compare the results that we obtain for the isoscalar coefficients a in Tab. IVII and those obtained 
in ref. [33|. The results of the fit 2+ of the Tab. 4 in ref. [13 are: 

2ci-c 3 = -1.63 ± 0.9 GoV" 1 

c 2 = -1.49 ± 0.67 GeV" 1 . (41) 



The agreement with the global fit from Tab. IVII is fair within errors. 
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TABLE VIII: Additional fits with D and different a nN 





D 3 a 


xi 


2.3 


D [fm] at 1110 MeV 


-0.011 


Ct-x-p 


-0.717 ±0.59 


a TV°n 


-0.823 ±0.26 


^TT + p 


-1.130 ±0.10 


2ci - c 3 [GeV- 1 ] 


-1.752 ±0.22 


c 2 [GeV" 1 ] 


-2.676 ±0.27 


/3 [MeV~ 2 ] 


0.002731 ± 1.8 ■ 10~ 4 


7 [10- 5 ■ m«] 


10.0 ±9.0 



A more direct comparison with the results of ref. |33j ] can be achieved by removing the damping factor j3 in the 
isoscalar term. We can only get good agreement with data up to about y/s = 1180 MeV. Restricting ourselves to 
energies below y/s = 1180 MeV, the fit to the data provides the parameters given in Tab. IVII as 'No damping'. One 
can sec that the C2 coefficient is in good agreement with eq. (|41|l and also the combination of 2c\ — C3 within errors. 
With this fit to the restricted data we find 

(6 , 61) = (-37 ± 37, -887 ± 43) • [10" V" 1 ]. (42) 

It is worth stressing this agreement since the "standard" values of the Cj coefficients used in chiral perturbative 
calculationSjWhere the A is not taken into account explicitly, are much larger in size than the coefficients found here 
or in ref. [33| and lead to the combinations of eq. I|41|) with opposite sign. Given the large amount of problems 
originated by the use of the "standard" (large) c, coefficients in chiral perturbative calculations, the models to fit 7^ 
cross sections including explicitly the A, as in ref. |3^|, and the "smaller" a coefficients found here and in ref. [33| . 
are highly recommendable. 

As for the 7 parameter corresponding to the irN — ► ttttN mechanism, we should expect on physical grounds quite 
a small contribution. Indeed, this is the case: The size of the 7rA^ — > 7rA^ term including two ttN — > ttttN vertices and 
the ttitN loop for the global fit 2 corresponds to about 5% of the tree level 7riV amplitude. 



6. Isospin breaking 

The isospin brea king in pion nucleon scattering has received much attention both phenomenologically and theo- 
retically [181, [13, Hfl lU I23, |23| . Our theoretical model uses isospin symmetry up to breaking effects from the use of 
different masses. We also rely for the fit upon some scattering data that has been analyzed assuming isospin symmetry 
|56j . A possible measure of isospin breaking can be given by the quantity D which describes the deviation from the 
triangle identity, 

D = /cex - (/+ - /-). (43) 

In ref. j^j, D is calculated by fitting s and p-wave of /+, /_, /cex which are the amplitudes of n + p, n~p, and 
ir~p — > ir°n. The authors fit a variety of potential models to the elastic channels, and predict /cex from that. Fitting 
in a second step /cex alone, they state a discrepancy between the prediction and the fit of /cex that results in a 
value of D = —0.012 ± 0.003 fm. Physical masses are included in the coupled channel approach of ref. [23], and their 
value for D contains effects of other sources of isospin violation than mass splitting. 

In the present work we have derived a microscopical model only for the s-wave interaction, also using physical 
masses in a coupled channels approach. Unfortunately, there is no partial wave analysis available that is free of 
isospin assumptions, so that we cannot evaluate D from eq. H43f) using only experimental s-wave amplitudes. For 
this reason we shall take advantage of the work done in [2(j • 

The global fit from Tab. | VII which only contains an isospin violation from mass splitting, produces D = —0.0066 fm, 
around half the value of ref. 20] . In order to give more freedom to the model to violate isospin symmetry from different 
sources than mass splitting, we now allow different subtraction constants in each channel, n~p, 7r°n, and ir + p. Then, 
we add extra data points taking the value of D from ref. |2(j at three energies covering the range of TJ a b = 30 — 50 
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MeV (as in Fig. 1 of ref. 20]). The x 2 stays practically the same compared to the global fit in Tab. I VII but the 
fit gives D = —0.011 fm at y/s — 1110 Mev, see 'D 3 a' in Tab. IVIIII The scattering amplitudes from this fit are 
practically the same as those shown in the Figs. and El for the global fit and we do not plot them again. 

As we can see in Tab. IVIIII the subtraction constant for the ir + p channel barely changes with respect to the global 
fit in Tab. I VII while those for the n~p and Tr°n channels are reduced in size by about 30 %. The values that we obtain 
in this fit for 2ci — C3 and C2 change little with respect to those quoted before and the values for (bo, b\) reveal a shift 
compared to the ones of the global fit in eq. H38JI of 

Sb a = 6 • 10~ 4 m-\ 5b l = 9 • 10~ 4 m" 1 . (44) 

The changes induced by the extra isospin breaking are rather small compared with the errors that we already have. 

The last value obtained for D is in agreement with ref. 113. A warning should be given about the solution found, 
since the isospin 1/2 amplitudes used in the fit imply isospin symmetry. However, the fact that half the amount of D 
that we obtain comes from the use of different masses without invoking isospin breaking from other sources, and that 
the threshold data, which have small error bars, and thus a large weight in our fit, do not imply isospin symmetry, 
makes us confident that the solution obtained accounts reasonably for isospin violation in the problem. 

V. CONCLUSIONS 

For a determination of the isoscalar and isovector scattering lengths of the ttN system, new calculations on the 
complex pion deuteron scattering length have been performed. The imaginary part of a^- d shows a very good 
agreement with experimental data. The dispersive part from absorption has been found to be compatible with zero. 
This, together with corrections from crossed diagrams and the A(1232) resonance, and with other corrections taken 
from the literature, leads to a substantial shift of the real part of a^-^ towards positive values. 

The unitary coupled channel approach of ref. |2j| has been tested for consistency at low energies. However, we 
have added an isoscalar term that can be matched to terms of higher orders of the chiral Lagrangians. These terms 
are known to play an important role at threshold. With this additional ingredient to the model, together with the 
ttN — > 2ttN channel, an acceptable global fit for the ttN amplitude up to intermediate energies has been obtained. 

One of the results of the present work concerns the values of the Cj parameters used in chiral perturbation theory 
at low energies. We find them compatible with values obtained from fits to data when the A is explicitly taken 
into consideration. On the other hand, we have addressed the isospin violation issue and found that our fit to data 
accounts for about half the isospin breaking only from mass splittings. The model has been extended to account for 
other sources of isospin breaking and then can match results of isospin breaking found in other works. We find that 
the effect of this breaking in the 60, b\ parameters is well within uncertainties from other sources. 

Attention to the sources of errors has been paid and we find larger values than in former studies. Altogether, we 
have here a new determination of the ttN scattering lengths and a new parametrization of the ttN s-wave amplitudes 
at low energies that can be used as input in studies of other elementary processes or as input to construct optical 
potentials from pionic atoms, where problems tied to the strength of the isoscalar part of the potential still remain. 
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APPENDIX: THE d-WAVE IN THE DEUTERON 

In the sections IIII BI and IIII CI the influence of the d-wave and the interference of d-wave and s-wave have been 
discussed for a number of diagrams, namely the absorption, and the 5th diagram of Fig. [S] We display some explicit 
formulas for the coupling of spin and angular momentum involved in these calculations. 

In momentum representation, the deuteron wave functions in eqns. (|ltj|) . H22fl . an d l|25|l are given by 

J= i p j 
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The index i stands for s or d-wave, and the parametrizations of the radial part ^(p), by means of the Cj, Dj, and 
nij, are taken from ref. |4l| ( n = H) f° r the CD-Bonn and ref. ^2] (n = 13) for the Paris wave function. The 
normalization is / dp p 2 (ty( s \p) 2 + ^^(p) 2 ) = 1 in both cases. 

The angular structure of the <i-wave is given by the angular momentum 1 = 2, l z = ji and the spin wave function \ 



C(2,l,l; M,f-^^2, M (0 P ,^)xi = 4-u , /^=^ooXs 



,S=1 



(A.2) 



The angular structure , also normalized to one, is preserved under Fourier transforms from coordinate space. The 
index v = —1, 0, 1 indicates at which 3rd component J z the total angular momentum is fixed. The C's in eq. (| A. 2|> 
are the Clebsch-Gordan coefficients that couple the nucleon spins to the angular momentum I = 2 of the d-state, in 
order to give a total angular momentum J of one. The nucleons are necessarily in a spin triplet with 5=1 for s and 
d-wave. 

Fixing in all calculations the third component of J at J z = v = 0, we obtain for the d — s interference of the 
diagrams of absorption in Fig. [3 



clcon 1.2 



q) {a* 



clcon 



1,2 -q') f(q + l)/r(q' + l 



1 



1-3 cos 2 



q+i 



2V27T 

3 

— =— cos 9 — r~r sin 9 — ~r 

2\/27r q+1 q+1 



q- q 



sm< 



q+i 



(q' x q) x - cos0^j(q' x q 



(A.3) 



where symmetries in the amplitude A in eq. (|22|l have been used. The a matrices act on the same nucleon. The 
second line corresponds to the first term in the decomposition erq' erq = q' • q + i (q' x q) • <r, and the third line to 
the term with crossed momenta. 

Another spin structure is given by the diagrams in Fig. [21 and the 5th diagram of Fig. [S] Here, the ttNN vertices 
are attached at different nucleons, and we obtain for the d — s interference for Fig. [3] 



(<?nucloon 1,2 ' q) (<?nucleon 2,1 ' q') fc > (q + 1) fo ' + 1) 
1 



f( s )/ 



2\/27r 
3 



1-3 cos 2 a 



q+i 



Izlz 



— =— cost/-— T smf-— v 
2V2n q+1 q+1 



q+i 



(A.4) 



The term {q x q' x + qyq'y ~ qzq' z ) m the second line of eq. (|A.4fl is also present in the s-wave — > s-wave transition and 
shows explicitly the 1/3-contribution of the diagrams in Fig. [3]to the ones of Fig. as has been derived in a different 
way in eqns. JTHJ, CHI), and PU|). 

For the 5th diagram of Fig. one obtains the angular structure by replacing the momentum components qi — > gr- 



and 



q+i 



of eq. 



U is 



in eq. ljA.4l) . where q' is the momentum of the disconnected pion (the q in the angles 6 } 
not changed). 

The angular structure of d-wave — > ti-wave transition for the diagrams is calculated in analogy to eqns. (|A.3(1 and 
(jA.4(l . but the resulting expressions are lengthier due to the occurring double sums from eq. (|A.2|I . 

In the calculations for absorption, the results have been numerically tested for choices of v equal to ±1 instead of 
0, and they stay the same, as it is required. 
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